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We present a direct measurement of the mass difference between top and antitop quarks (Am) in lepton+jets 
tt final states using tiie "matrix element" method. The purity of the lepton+jets sample is enhanced for tt events 
by identifying at least one of the j et as originating from a b quark. The analyzed data correspond to 3 . 6 fb~ ^ of 
pp collisions at ^/s = 1.96 TeV acquired by DO in Run II of the Fermilab Tevatron Collider. The combination 
of the e+jets and /i+jets channels yields Am = 0.8 ±1.8 (stat) it 0.5 (syst) GeV, which is in agreement with 
the standard model expectation of no mass difference. 

PACS numbers: 14.65.Ha 



I. INTRODUCTION 

The standard model (SM) is a local gauge-invariant quan- 
tum field theory (QFT), with invariance under charge, parity, 
and time reversal (CPT) providing one of its most fundamen- 
tal principles [1-4], which also constrains the SM [5]. In 
fact, any Lorentz-invariant local QFT must conserve CPT [6]. 
A difference in the mass of a particle and its antiparticle 
would constitute a violation of CPT invariance. This issue has 
been tested extensively for many elementary particles of the 
SM [7]. Quarks, however, carry color charge, and therefore 
are not observed directly, but must first hadronize via quan- 
tum chromodynamic (QCD) processes into jets of colorless 
particles. These hadronization products reflect properties of 
the initially produced quarks, such as their masses, electric 
charges, and spin states. Except for the top quark, the time 
scale for hadronization of quarks is orders of magnitude less 
than for electroweak decay, thereby favoring the formation of 
QCD-bound hadronic states before decay. This introduces a 
significant dependence of the mass of a quark on the model of 
QCD binding and evolution. In contrast to other quarks, no 
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bound states are formed before decay of produced top quarks, 
thereby providing a unique opportunity to measure directly 
the mass difference between a quark and its antiquark [8]. 

In proton-antiproton collisions at the Fermilab Tevatron 
Collider, top quarks are produced in tt pairs via the strong 
interaction, or singly via the electroweak interaction. In the 
SM, the top quark decays almost exclusively into a W boson 
and a b quark. The topology of a tt event is therefore deter- 
mined by the subsequent decays of the W bosons. The world's 
most precise top quark mass measurements are performed in 
the lepton-|-jets (£-|-jets) channels, which are characterized by 
the presence of one isolated energetic electron or muon from 
one W ^ £v decay, an imbalance in transverse momentum 
relative to the beam axis from the escaping neutrino, and four 
or more jets from the evolution of the two b quarks and the 
two quarks from the second W qq' decay. 

The top quark was discovered [9, 10] in proton-antiproton 
collision data at a center of mass energy of = 1.8 TeV in 
Run I of the Tevatron. After an upgrade to a higher center 
of mass energy of y/s = 1 .96 TeV and higher luminosities. 
Run II of the Tevatron commenced in 2001. Since then, a 
large sample of ti events has been collected, yielding precision 
measurements of various SM parameters such as the mass of 
the top quark, which has been determined to an accuracy of 
about 0.6% or mtop = j{m,+mr) = 173.3 ± 1.1 GeV [11], 
where mt {mi) is the mass of the top (antitop) quark. 

The DO Collaboration published the first measurement of 
the top-antitop quark mass difference, I^m = nit — mj, using 
1 fb ' of Run II integrated luminosity [12]. Our new measure- 
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ment, presented here, employs the same matrix element (ME) 
technique [13, 14], suggested initially by Kondo et al. [15- 
17], and developed to its current form by DO [18]. Our previ- 
ous study measured a mass difference 

Am = 3.8 ± 3.4 (stat.) ± 1.2 (syst.) GeV. 

Recently, CDF has also measured Am [19] based on 5.6 fb~' 
of Run II data, using a template technique, and found 

Am = -3.3 ± 1.4 (stat.) ± 1.0 (syst.) GeV. 

In this paper, we extend our first measurement of Am us- 
ing an additional 2.6 fb^ ' of Run II integrated luminosity, and 
combining our two results. We also re-examine the uncer- 
tainties from the modeling of signal processes and of the re- 
sponse of the detector. Most important is a possible presence 
of asymmetries in the calorimeter response to b and fo-quark 
jets, which we re-evaluate using a purely data-driven method. 
We also consider for the first time a bias from asymmetries in 
response to c and c-quark jets. 

This paper is arranged as follows: after a brief descripton of 
the DO detector in Sec. II, we review the event selection and 
reconstruction in Sec. III. In Sec. IV, we define the samples of 
Monte Carlo (MC) events used in the analysis. The extraction 
of the top-antitop quark mass difference using the ME tech- 
nique is then briefly reviewed in Sec. V. The calibration of this 
technique, based on MC events, and the measurement of the 
mass difference in 2.6 fb~^ of Run II integrated luminosity are 
presented in Sec. VI. The evaluation of systematic uncertain- 
ties and cross checks are discussed in Sec. VII and VII C, re- 
spectively. Finally, the combination of the measurements for 
the 2.6 fb~^ and 1 fb~^ data samples is presented in Sec. VIII. 

II. THE DO DETECTOR 

The DO detector has a central-tracking system, calorime- 
try, and a muon system. The central-tracking system con- 
sists of a silicon microstrip tracker (SMT) and a central fiber 
tracker (CFT), both located within a 1.9T superconducting 
solenoidal magnet [20-22], with designs optinnized for track- 
ing and vertexing at pseudorapidities j?] | < 3 [23]. The SMT 
can reconstruct the pp interaction vertex (PV) with a precision 
of about 40 nm in the plane transverse to the beam direction 
and determine the impact parameter of any track relative to 
the PV [24] with a precision between 20 and 50 jUm, depend- 
ing on the number of hits in the SMT. These are the key ele- 
ments to lifetime-based Z^-quark jet tagging. The liquid-argon 
and uranium sampling calorimeter has a central section (CC) 
covering pseudorapidities |tj | < 1.1 and two end calorimeters 
(EC) that extend coverage to |t7| « 4.2, with all three housed 
in separate cryostats [20, 25]. Central and forward preshower 
detectors are positioned just before the CC and EC. An outer 
muon system, at 1 77 1 < 2, consists of a layer of tracking detec- 
tors and scintillation trigger counters in front of 1.8 T toroids, 
followed by two similar layers after the toroids [26]. The lu- 
minosity is calculated from the rate of pp inelastic collisions 
measured with plastic scintillator arrays, which are located in 



front of the EC cryostats. The trigger and data acquisition 
systems are designed to accommodate the high instantaneous 
luminosities of Run II [27]. 



ni. EVENT SELECTION 

In this new measurement of Am, we analyze data corre- 
sponding to an integrated luminosity of about 2.6fb~^ for both 
the e-|-jets and ;U-|-jets channels. 

Candidate tt events are required to pass an isolated ener- 
getic lepton trigger or a lepton+jet(s) trigger. These events 
are enriched in tt content by requiring exactly four jets re- 
constructed using the Run II cone algorithm [28] with cone 
radius AR = ^ (Ar))^ + (A^)^ = 0.5, transverse momenta 
Pt > 20 GeV, and pseudorapidities |t7| < 2.5. The jet of 
highest transverse momentum in a given event must have 
Pt > 40 GeV. Furthermore, we require exactly one isolated 
electron with pj > 20 GeV and |?7| < 1.1, or exactly one iso- 
lated muon with pr > 20 GeV and |tj| < 2.0. The leptons 
must originate within 1 cm of the PV in the coordinate along 
the beamline. Events containing an additional isolated lep- 
ton (either e or 11) with pr > 15 GeV are rejected. Lepton 
isolation criteria are based on calorimetric and tracking infor- 
mation along with object identification criteria, as described 
in Ref. [29]. The positively (negatively) charged leptons are 
used to tag the top (antitop) quark in a given event. To reduce 
instrumental effects that can cause charge-dependent asym- 
metries in the lepton momentum scale, the polarity of the 
solenoidal magnetic field is routinely reversed, splitting the 
total data into two samples of approximately equal size. The 
PV must have at least three associated tracks and lie within the 
fiducial region of the SMT. At least one neutrino is expected 
in the £+jets final state; hence, an imbalance in transverse 
momentum (defined as the opposite of the vector sum of the 
transverse energies in each calorimeter cell, corrected for the 
energy carried by identified muons and energy added or sub- 
tracted due to the jet energy scale calibration described be- 
low) of > 20 GeV (25 GeV) must be present in the e-|-jets 
(jU+jets) charmel. These kinematic selections are sunomarized 
in Table 1. 

To reduce the contribution of multijet production (MJ) in 

the e+jets channel, A(p(e,i^T) > 2.2- j^r x 0.045 GeV"^ 
is required for the azimuthal difference A(j){e,^T) = \(l>e — 
<j)^j.\ between the electron and the direction of ^r- Like- 
wise, A<j){p.,i^T) > 2.1 — j^T 'X 0.035 GeV"' is required in 
the /i+jets channel. Jets from quarks are identified by a 
neural-network-based ^-tagging algorithm [30], which com- 
bines variables that characterize properties of secondary ver- 
tices and tracks within the jet that have large impact parame- 
ters relative to the PV. Typically, its efficiency for ^-quark jets 
is about 65%, while the probability for misidentifying u, d, 
s-quark and gluon jets as b jets is about 3%. To increase tt 
purity, and to reduce the number of combinatoric possibilities 
for assigning jets to tt decay products, we require at least one 
fo-tagged jet to be present in the events used to measure Am. 

After all acceptance requirements, a data sample of 312 
(303) events is selected in the e-|-jets (;U-|-jets) channel. As 
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TABLE 1 : A summary of kinematic event selections applied. 



Exactly 1 charged lepton 



pr>20GeV |t]| < 1.1 (e) 

PT >2QG€W \^\<2.Q{^i) 

Exactly 4 jets pr > 20 GeV | t] | < 2.5 

Jet of highest pj Pt > 40 GeV | t] | < 2.5 

Imbalance in transverse momentum ^ 

> 25 GeV (ju+jets) 



discussed above, each of those samples is split according to 
lepton charge. In the e+jets channel, 174 (138) events have 
a positive (negative) lepton in the final state. Likewise, the 
/i+jets sample is split to subsets of 145 and 158 events. 



We check our modified PYTHIA version against the orig- 
inal by comparing large samples of simulated tt events for 
{mt,mt) = (170 GeV, 170 GeV), at both the parton and re- 
construction levels, and find full consistency. 

The tt samples are generated at fourteen combinations of 
top and antitop quark masses {mt,mf), which form a grid 
spaced at 5 GeV intervals between (165 GeV, 165 GeV) and 
(180 GeV, 180 GeV), excluding the two extreme points at 
(165 GeV, 180 GeV) and (180 GeV, 165 GeV). The four 
points with mt = are generated with the standard pythia, 
whereas all others use our modified version of the generator. 



B. Monte Carlo and other simulations of baclcground 



rV. MONTE CARLO SIMULATION 

Large samples of simulated MC events are used to de- 
termine the resolution of the detector and to calibrate the 
Am measurement as well as the statistical sensitivity of the 
method. After simulation of the hard scattering part of the in- 
teraction and parton shower corrections, MC events are passed 
through a detailed detector simulation based on GEANT [31], 
overlaid with data collected from a random subsample of 
beam crossings to model the effects of noise and multiple in- 
teractions, and reconstructed using the same algorithms that 
are used for data. Although the fraction of signal events, /, is 
fitted in the analysis, we also cross check that the entire data 
sample is described adequately by the simulations. 



A. Monte Carlo samples for signal 

Simulated tt events with different mt and mj are required to 
calibrate the Am measurement. We use the PYTHIA genera- 
tor [32], version 6.413, to model the ti signal. This generator 
models the Breit-Wigner shape of the invariant mass distribu- 
tion of t and t quarks, whose correct description is important 
for the Am measurement. 

In the standard pythia, it is not possible to generate tt 
events with different masses mt and mf. Therefore, we modify 
the PYTHIA program to provide signal events with mt ^ mf. 
In applying these modifications, we adjust the description of 
all quantities that depend on the two masses, for example, the 
respective decay widths F, and F,. Technical details of this 
implementation can be found in Appendix 1. 

We generate ti events using the CTEQ6L1 parton distribu- 
tion function set (PDF) [33] at the momentum transfer scale 

transverse momentum for the hard scattering process, and f, 
is the four-momentum of the incoming parton i. For mt = mj, 
the expression used for is identical to that in the standard 
PYTHIA. All other steps in the event simulation process aside 
from the generation of the hard- scattering process, e.g., the 
modeling of the detector response, are unchanged from the 
standard pythia. 



The dominant background to ti decays into £-|-jets final 
states is from the electroweak production of a W boson in 
association with jets from gluon radiation. We simulate the 
hard scattering part of this process using the ALPGEN MC pro- 
gram [34], which is capable of simulating up to five additional 
particles in the final state at leading order (LO) in «(. ALPGEN 
is coupled to pythia, which is used to model the hadroniza- 
tion of the partons and the evolution of the shower. The 
MLM matching scheme is applied to avoid double-counting 
of partonic event configurations |35]. The W+jets contribu- 
tion is divided into two categories according to parton flavor: 
(/) W +bb+ jets and W+cc +jets and (ii) all other contribu- 
tions, where "jets" generically denotes jets from M,c/,s-quarks 
and gluons. The second category also includes the W+c+jets 
final states. While the individual processes are generated with 
ALPGEN, the relative contributions of the two categories are 
determined using next-to-LO (NLO) calculations, with next- 
to-leading logarithmic (NLL) corrections based on the MCFM 
MC generator [36]. These NLO corrections increase the LO 
cross section of category (/) by a factor of k = I Al ±0.22, 
while A: = 1 is used for category (//). The resulting combined 
W+jets background contribution is then determined from a fit 
to data and predictions for other signal and background con- 
tributions, as described in Sec. V. Thus, the NLO fc-factors 
only change the relative balance between (;) and (//). 

Additional background contributions arise from WW, WZ, 
ZZ, single top quark electroweak production, Z tt, and 
Z ^ ee (Z ^ jiji) production in the e+jets (/i+jets) chan- 
nel. The predictions for these backgrounds are taken from 
MC simulations, and, with the exception of single top quark 
electroweak production, their production cross sections are 
normalized to NLO+NLL calculations with MCFM. Diboson 
processes are simulated with pythia. The hard- scattering 
part of single top quark production is simulated with COM- 
pHEP [37], while alpgen is used for Z-|-jets boson produc- 
tion. For both backgrounds, PYTHIA is employed to model 
hadronization and shower evolution. The CTEQ6L1 PDFs 
and the DO Tune A underlying event model [38] are used in 
the generation of all MC samples. 

Events from MJ production can pass our selection crite- 
ria, which typically happens when a jets mimics an electron, 
or a muon that arises from a semileptonic decay of a /? or 
c quark appears to be isolated. The kinematic distributions 
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of the MJ background are modeled using events in data that 
fail only the electron identification (muon isolation) criteria, 
but pass loosened versions of these criteria defined in [40]. 
The absolute contribution of this background to each of the 
channels is estimated using the method described in Ref. [40]. 
This method uses the absolute numbers of events with prompt 
leptons A^ioos^ and events from MJ production Nf^^^^ in the 
sample with loosened lepton identification criteria, and relates 
them to the absolute contributions to the sample with standard 
lepton identification criteria via = e''+^N[l+^ + e'^^N^^^^. 
Here, e"+'^ and e'^^ represent the efficiency of events which 
pass the loosened lepton identification criteria to also pass the 
standard identification criteria, and are measured in control 
regions dominated by prompt leptons and MJ events, respec- 
tively. 



C. Event yields 

We split the selected £+jets events into subsamples accord- 
ing to lepton flavor (e or jj.), jet multiplicity, and the number of 
/7-tagged jets in the event to verify an adequate description of 
the data with our signal and background model. In general, we 
observe good agreement between data and simulations, and 
systematic uncertainties on the final result explicitly account 
for moderate agreement observed in some kinematic distribu- 
tions (cf. Sec. VII). 

The numbers of events surviving the final stage of selection 
with at least one b-tag are summarized in Table 2. Here, for 
ease of comparison, the contributions from tt events are scaled 
to 7.45+0 7 pb, section including NNLO ap- 

proximations [41]. The total W+jets cross section is adjusted 
to bring the absolute yield from our signal and background 
model into agreement with the number of events selected in 
data before applying ^-jet identification criteria. The distri- 
butions in the transverse mass of the W boson, [42], and 
in 7- are shown in Fig. 1 for data with at least one b-tag, to- 
gether with the predictions from our signal and background 
models. 



V. GENERAL DESCRIPTION OF THE METHOD 

In this section, we describe the measurement of Am using 
the ME method. The procedure is similar to the one used 
in Ref. [13, 43] to measure the average top quark mass mtop, 
but instead of simultaneously determining mtop and the jet en- 
ergy scale (JES), here we measure directly the masses of the 
top and antitop quarks, m, and m,-, which provides Am and 
mtop. We review the ME approach in Sec. V A, the calcula- 
tion of signal and background event probabilities in Sees. V B 
and V C, respectively, as well as the parametrization of the 
detector response and the use of /j-tagging information in 
Sec. VD. 



TABLE 2: Numbers of events selected in data, compared to yield 
predictions for individual processes using simulations, in the e+jets 
and ^+jets channels with exactly 4 jets and at least one Z?-tagged jet, 
split according to i-tag multiplicity. Uncertainties are purely statis- 
tical. See text for details. 





1 /7-tag 


^ 1 /?-t3gS 


e+jets 






tt 




VI. ± Z.J 


W +jets 


iy.y ± l.z 


4. / ± U.J 


iVlJ 


1^ c 1 1 
Zj.J ± Z. 1 


J. / ±: i.U 


Z+jets 


7.6 ± 0.7 


0.9 ±0.1 


Other 


6.6 ± 0.4 


1.9 ±0.1 


Total 


216.7 ± 3.9 


105.1 ± 2.7 


Observed 


223 


89 


/i+jets 






tt 


105.9 ± 2.4 


70.9 ± 2.0 


ly+jets 


59.9 ± 1.8 


7.2 ± 0.5 


MJ 


5.2 ± 0.9 


2.0 ± 0.6 


Z+jets 


5.3 ± 0.5 


1.2 ±0.2 


Other 


5.0 ±0.3 


1.3 ±0.1 


Total 


181.3 ± 3.2 


82.6 ± 2.2 


Observed 


191 


112 




FIG. 1 : The transverse mass of the W boson for events with at 
least one ii-tag is shown for the (a) e+jets and (b) ^u+jets channels. 
Similarly, is shown for the (c) e+jets and (d) /J+jets channels. 
The statistical uncertainties on the prediction from the tt signal and 
background models are indicated by the hatched area. 



A. Probability densities for events 

To optimize the use of kinematic and topological informa- 
tion, each event is assigned a probability Pgvi to observe it 
as a function of the assumed top and antitop quark masses: 
^evt = ^evt('«r,'«f)- The individual probabilities for all events 
in a given sample are combined to form a likelihood, from 
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which the Am and Wtop parameters are extracted. Simphfy- 
ing assumptions are made in the expression of the likelihood 
about, e.g., detector response or the sample composition, are 
made to render the problem numerically solvable. It is there- 
fore necessary to calibrate the method using fully simulated 
MC events, as detailed in Sec. VI B. Systematic uncertainties 
are estimated to account for possible effects of these assump- 
tions on the extracted value of Am. 

Assuming that the signal and background physics processes 
do not interfere, the contribution to the overall probability 
from a single event can be formulated as 



2qiq2S 



dOe, (3) 



Penix;m,,mt,f) 



f ■Psig{x;mt,mt) 
+ (l-/)-/lkgW}, 



(1) 



where x denotes the set of measured kinematic variables for 
the event observed in the detector, / is the fraction of sig- 
nal events in the sample, A{x) reflects the detector acceptance 
and efficiencies for a given x, and Psig and /^kg are the prob- 
abilities for the event to arise from tt or W+jets production, 
respectively. The production of W bosons in association with 
jets is the dominant background, and we neglect all other con- 
tributions to Pbkg- Kinematically similar contributions from 
other background processes like MJ production are accounted 
for in the analysis implicitly (cf. Sec. VII). 

Both signal and background probabilities depend on the 
JES, which is defined as the ratio of the calibrated energy of a 
jet over its uncahbrated energy. The standard calibration of jet 
energies accounts for the energy response of the calorimeters, 
the energy that crosses the cone boundary due to the transverse 
shower size, and the additional energy from pileup of events 
and from multiple pp interactions in a single beam crossing. 
Although the Am observable is not expected to show a strong 
dependence on JES by construction, we apply an additional 
absolute calibration to the JES using a matrix element which 
is a function of ;7ttop and JES from Refs. [13, 43]. The poten- 
tial systematic bias on Am from the uncertainty on the absolute 
value of the JES is estimated in Sec. VII. 

To extract the masses mt and mf from a set of n selected 
events, with sets of measured kinematic quantities xi,...,x„,a 
likelihood function is defined from the individual event prob- 
abilities according to Eq. (1): 

n 

L{xi,...,x„;m!,mj,f) = Y[Pevt{xi; mt,mf,f) . (2) 

i=l 

For every assumed {mi,mf) pair, we first determine the value 
of / = /''^^' that maximizes this likelihood. 



B. Calculation of signal probability Pgig 

The probability density for the signal to yield a given set of 
partonic final state four-momenta y in pp coUisions is propor- 
tional to the differential cross section da for tt production: 

do{pp^ti^y;mt,mt)= / ^ d^id^2/(^i)/(^2) 

quark 
11 '92 flavors 



where ^ denotes the matrix element for the qq ^ tt ^ 
b{lv)b{qq') process, s is the square of the center-of-mass en- 
ergy, qi is the momentum fraction of the colliding parton / (as- 
sumed to be massless), and d<t>6 is an infinitesimal element of 
six-body phase space. The f{qi ) denote the probability densi- 
ties for finding a parton of given fiavor and momentum frac- 
tion qi in the proton or antiproton, and the sum runs over all 
possible flavor configurations of the colliding quark and an- 
tiquark. In our definition of and therefore the tt signal 
probability, only quark-antiquark annihilation at LO is taken 
into account; in this sense, Eq. (3) does not represent the full 
differential cross section for tt production in pp collisions. 
Effects from gluon-gluon and quark-gluon induced ti produc- 
tion are accounted for in the calibration procedure described 
in Sec. VI B. We further test for an effect on Am from higher- 
order corrections in Sec. VII C. 

The differential cross section for observing a tt event with 
a set of kinematic quantities x measured in the detector can be 
written as 

Aa{pp ^ tt ^ x; mt , mj, fcjEs ) 

= A{x) dy do {pp^tt^y;mt, mt)W{x, y; kms ) , (4) 

Jy 

where finite detector resolution and offline selections are taken 
explicitly into account through the convolution over a transfer 
function W{x,y;kj£s) that defines the probability for a par- 
tonic final state y to appear as x in the detector given an abso- 
lute JES correction ^jes- 

With the above defintions, the differential probability to ob- 
serve a tt event with a set of kinematic quantities x measured 
in the detector is given by 

D / , N da{pp^tt^x;mt,mf,k3Es) 
Psig[x; m,, nit, k JES) = . _ , ; (5) 

where Cobs is the cross section for observing tt events in the 
detector for the specific ME ^ defined in Eq. (3): 

O'obs (PP^ tt; mt , mf, fcjES ) 

= / dxdy da{pp ti y;mt,mf)W{x,y;k3Es)Mx) 



= dy da{pp -^tt^ y;mt,mt) dxW{x,y;kjEs)A{x) . 

Jy Jx 

The normalization factor aobs is calculated using MC integra- 
tion techniques: 

OobsiPP ^ tt;mt,mt,k]Es) - atot{mt,mt)x{A\mt,mt), (6) 
where 

atot {mt , mj) = dy do{pp ^ tt ^ y;mt, m, ) , (7) 

Jy 

and 



(A|m;,mf) = — £ 

iVc 



(0. 



(8) 



gen acc 
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To calculate the {A\mt,mf) term, events are generated ac- 
cording to da{pp tt;mt,mf) using PYTHIA and passed 
through the full simulation of the detector Here, A^gen is the 
total number of generated events, co are the MC event weights 
that account for trigger and identification efficiencies, and the 
sum runs over all accepted events. 

The formulae used to calculate the total cross section CTtot 
and the matrix element are described below in Sees. V B 1 
and V B 2. In all other respects, the calculation of the sig- 
nal probability proceeds identically to that in Refs. [13, 43], 
with the following exceptions: (/) CTEQ6L1 PDFs are used 
throughout, and (//) the event probabilities are calculated on a 
grid in m, and nif spaced at 1 GeV intervals along each axis. 
As described in Sec. VI A, a transformation of variables to Am 
and mtop is performed when defining the likelihood. 



1. Calculation of the total cross section Otot 

Without the assumption of equal top and antitop quark 
masses, the total LO cross section for the qq tt process in 
the center of mass frame is given by 



21 si 



(9) 



where Et (Ef) are the energies of the top and antitop quark, 
and p is the three-momentum of the top quark. This reduces 
to the familiar form for m, = mf. 



9s ^\ 3 



where )3 = \pt\/Et — \pf\/Ef represents the velocity of the t 
(or f) quark in the qq rest frame. 

Integrating Eq. (9) over all incoming qq momenta and using 
the appropriate PDF yields Oioi{pp — > tt; mt,mf), as defined 
for any values of nit and m,- in Eq. (7). Figure 2 displays the 
dependence of atot on Am for a given mtop. The correspond- 
ing average acceptance term (A|m(,m,-), as defined in the same 
equation, is shown in Fig. (3) for the e+jets and /x+jets chan- 
nels. 



2. Calculation of the matrix element ^ 

The LO matrix element for the qq tt process we use in 
our analysis is 



l±FF- 
9 s 



:{(£,- \p,\c,i, f + {Er+ \pj\c^tf +2m,mj] .(10) 



The form factors FF are identical to those given in Eqs. (24) 
and (25) of Ref. [13]. For the special case of m, = m,, the 
expression in Eq. (10) reduces to 



\^\^^^-±FF-{2~fih],) 



Q. 

^ 11 

o 

'^^'lO 

g 

8 
7 
6 
5 
4 
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FIG. 2: The total pp tt production cross section Otot defined in 
Eq. (7) as a function of Am and nitop- Each line shows Otot as a 
function of Am for a given value of mtop displayed above the curve. 
The range from 152 GeV to 188 GeV is shown in 6 GeV increments, 
the broken line corresponds to 170 GeV. 




20 40 
Am (GeV) 



20 40 
Am (GeV) 



FIG. 3: The dependence of the overall average acceptance {A\m,,mf) 
on Am and mtop, as defined in Eq. (8), for the (a) e+jets and (b) /x + 
jets signal MC samples. Each line shows {A\mt,mj) as a function of 
Am for a given value of mtop displayed above the curve. The range 
from 152 GeV to 188 GeV is shown in 6 GeV increments, the broken 
lines correspond to 170 GeV. 



which is identical to Refs. [13, 44], where s^t is the sine of 
the angle between the incoming parton and the outgoing top 
quark in the qq rest frame. 



C. Calculation of the background probability i^kg 

The expression for the background probability Pbkg is sim- 
ilar to that for fjig in Eq. (5), except that the ME ^vK+jets is 
for W +jets production, and all jets are assumed to be light 
quark or gluon jets. Clearly, ^vv+jets does not depend on m, 
or m,-, and P^kg is therefore independent of either. We use a 
LO parameterization of ./^ from the VECBOS [45] program. 
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More details on the calculation of the background probabiUty 
can be found in Ref. [13]. 

D. Description of detector response 

The transfer function W{x,y,kjEs), which relates the set of 
variables x characterizing the reconstructed final-state objects 
to their partonic quantities y, is crucial for the calculation 
of the signal probability according to Eq. (5), and the cor- 
responding expression for /^kg. A full simulation of the de- 
tector would not be feasible for calculating event probabilities 
because of the overwhelming requirements for computing re- 
sources. Therefore, we parametrize the detector response and 
resolution through a transfer function. 

In constructing the transfer function, we assume that the 
functions for individual final-state particles are not correlated. 
We therefore factorize the transfer function into contributions 
from each measured final-state object used in calculating P^ig, 
that is the isolated lepton and four jets. The poorly measured 
imbalance in transverse momentum ^t, and consequently the 
ttansverse momentum of the neutrino, is not used in defining 
event probabilities. We assume that the directions of e, pL, 
and jets in (rj , ^) space are well-measured, and therefore de- 
fine the transfer functions for these quantities as 5 functions: 
5^(77,^) = 5{r]y — r]x)5{^y — (px). This reduces the number of 
integrations over the 6-particle phase space dOg by 5 x 2 = 10 
dimensions. The magnitudes of particle momenta \p\ display 
sigruficant variations in resolution for leptons and jets and are 
therefore parameterized by their corresponding resolutions. 

There is an inherent ambiguity in assigning jets recon- 
structed in the detector to specific partons from ti decay. Con- 
sequently, all 24 permutations of jet-quark assignments are 
considered in the analysis. The inclusion of /j-tagging infor- 
mation provides improved identification of the correct per- 
mutation. This additional information enters the probability 
calculation through a weight Wi on a given permutation ; of 
jet-parton assignments. The w, are larger for those permu- 
tations that assign the Z?- tagged jets to b quarks and untagged 
jets to light quarks. The sum of weights is normalized to unity: 

Based on the above, we define the transfer function as 
W{x,y; kjEs) = We{E„Ey)5}{r],(l)) 

24 r 4 ^ 

^ L^' |n 52.(77, 0)Wjet(£;,£^^;/=jEs)|, (11) 

where £ denotes the lepton flavor, with a term We describing 
the energy resolution for electrons and the resolution in the 
transverse momentum for muons. Similarly, Wjet describes the 
energy resolution for jets. The sum in i is taken over the 24 
possible permutations of assigning jets to quarks in a given 
event. More details on Wg and Wjct can be found in Ref. [43]. 

The weight w,- for a given permutation / is defined by 
a product of individual weights w- for each jet j. For b- 
tagged jets, wj is equal to the per-jet tagging efficiency 
£tag{(^k', Ej, 7]'), where labels the three possible parton- 
flavor assignments of the jet: (?) b quark, (ii) c quark, and 



(//;) light {u,d,s) quark or gluon. For untagged jets, the 
factors are equal to 1 — etag(o!i; Ej,r]^). 

Because the contributions to IV+jets are parameterized by 
■^w+jets without regard to heavy-flavor content, the weights 
Wi for each permutation in the background probability are all 
set equal. 



VI. MEASUREMENT OF THE TOP-ANTITOP QUARK 
MASS DIFFERENCE 

A. Fit to the top-antitop quarlc mass difference 

For the set of selected events, the likelihood L{mt,mt) is 
calculated from Eq. (2) (Sec. V A). The signal fraction f^^^^ 
that maximizes the likelihood is determined at each {nit, mi) 
point for grid spacings of 1 GeV. Subsequently, a trans- 
formation is made to the more appropriate set of variables 
(Am,mtop): 

L(xi, ...,x„;AOT,OTtop) 

= L[xi,...,x„; AOT,OTtop,/''^"(AOT,mtop)]. (12) 

To obtain the best estimate of Am in data, the two- 
dimensional likelihood in Eq. (12) is projected onto the Am 
axis, and the mean value {h.m), that maximizes it, as well as 
the uncertainty on (Am) are calculated. This procedure 
accounts for any correlations between Am and mtop. As a con- 
sistency check, we simultaneously extract the average mass 
»Jtop by exchanging Am ^ m^p above. 



B. Calibration of the method 

We caUbrate the ME method by performing 1000 MC 
pseudo-experiments at each input point {nit, mi). These are 
used to correlate the fitted parameters with their true input 
values and to assure the correctness of the estimated un- 
certainties. Each pseudo-experiment is formed by drawing 
A^sig signal and A'bkg background events from a large pool of 
fully simulated ti and W-l-jets MC events. We assume that 
IV+jets events also represent the kinematic distributions ex- 
pected from MJ production and other background processes 
with smaller contributions, and evaluate a systematic uncer- 
tainty from this assumption. Events are drawn randomly and 
can be used more than once, and an "oversampling" correc- 
tion [46] is applied. The size of each pseudo-experiment, 
^ = ^sig + ^bkg. is fixed by the total number of events ob- 
served in the data, i.e., A'^ = 312 (303) events for the e-|-jets 
(/i+jets) channel. The fraction of signal events is allowed to 
fluctuate relative to the signal fraction / determined from data 
(Sec. VI B 1), assuming binomial statistics. The same W-|-jets 
background sample is used to form pseudo-experiments for 
each (m;,mf) mass point. 
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1. Determining the signal fraction in data 

The signal fraction / is determined independently for the 
e+jets and /i+jets channels directly from the selected data 
sample. The likelihood depends explicitly on three parame- 
ters: A;w, mtop, and /, as defined in Eq. (12). The uncalibrated 
signal fraction y""'^''' is calculated in data as an average of f°^^^ 
determined at each point in the [nit , m, ) grid and weighted by 
the value of the likelihood at that point. To calibrate y""'^^', 
we form 1000 pseudo-experiments for each input signal frac- 
tion /"^"^ in the interval [0, 1] in increments of 0. 1, and extract 
one, following the same procedure as in data. 
Signal MC events with m, = m,- = 172.5 GeV are used for this 
calibration. A linear dependence is observed between Z*^"" 
and where is the average of y""'^^' values extracted 
in 1000 pseudo-experiments for a given We use the re- 
sults of a linear fit of to f^"^^ to calibrate the fraction of 
signal events in data. The results are summarized in Table 3. 
Possible systematic biases on the measured value of Am from 
the uncertainty on / are discussed in Sec. Vll. 

TABLE 3: Signal fractions determined from data for the assumption 
that rut =mj = 172.5 GeV. The uncertainties are statistical only. 

Channel Measured signal fraction 
e+jets 0.71 ± 0.05 

/x+jets 0.75 ± 0.04 



TABLE 4: Fit parameters for the calibration of Am and mtop, defined 
by Eq. (13), and average pull- widths {(On) for pulls in Am and mtop, 
defined in Eq. (14). 





Channel 


<?0 (GeV) 


^1 


{(Ox) 


Am 


e^ 


-jets 


0.28 ± 0.14 


1.10 ±0.02 


1.25 ± 0.01 




fjets 


-0.08 ±0.13 


0.99 ± 0.02 


1.22 ±0.01 


mtop 


e^ 


-jets 


0.53 ± 0.08 


0.99 ± 0.02 


1.17 ± 0.01 




fjets 


0.24 ± 0.07 


1.02 ± 0.02 


1.16 ±0.01 



the 14 pull widths in each channel to constant offsets and are 
summarized in Table 4. We calibrate the estimated uncertainty 
according to = (a),r^„) x 5a,„. 



3. Calibration o/mtop 

Results from an analogous calibration of mtop are displayed 
in Fig. 5 (a) and (b) for the e+jets and ^+jets channel, respec- 
tively. The distributions in pull widths are given in parts (c) 
and (d) of the same figure. The corresponding fit parameters 
and average pull widths are also summarized in Table 4. 



2. Calibration of Am 

The dependence of the extracted Am on the generated Am 
is determined from the extracted values Am^^^' [mf^mf], again 
obtained from averaging (Am) over 1000 pseudo-experiments 
for each (m, , m,- ) combination. The resulting distribution and 
fit to the 14 (m,,m,-) points is shown in Fig. 4 (a) and (b) for 
the e+jets and /i+jets channels, respectively. This provides 
the calibration of the extracted Am value: 



Am'' 



^Am ^ ^ Am . ^^gen 



(13) 



The fit parameters ^ A'" are summarized in Table 4. 

For an unbiased estimate of Am and of the uncertainty ^Am 
on the measured (Am) value, the distribution of the pulls 
should be described by a Gaussian function with a standard 
deviation (SD) of unity, and centered at zero. A SD of the 
pulls larger than unity would indicate an underestimation of 
5a;,7, which could be caused by the simplifying assumptions 
of the ME technique discussed in Sec. V. For a given pseudo- 
experiment at (m,,m,-), we define the pull in Am as 



TTa,,, - 



{Am) - Am'"'"' (m, , mj) 



(14) 



"Am 



The pull widths (O-ji^^^, defined by the SD in Gaussian fits to the 
pull distributions, are also shown for all 14 (m,,m,-) points in 
Fig. 4 (c) and (d) for the e+jets and jU+jets channels, respec- 
tively. The average pull widths {co,i^J are taken from fits of 
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FIG. 4: The calibration of the extracted Am value as a func- 
tion of generated Am is shown for the (a) e+jets and (b) fl + 
jets channels. The points are fitted to a linear function. Each 
point represents a set of 1000 pseudo-experiments for one of 
the fourteen {mi,mj) combinations. The circle, square, trian- 
gle, rhombus, cross, star, and "x" symbols stand for mtop = 
165, 167.5, 170, 172.5, 175, 177.5, and 180 GeV, respectively. Sim- 
ilarly, the pull widths, as defined in the text, are given for the 
(c) e+jets and (d) /i+jets channels. 
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FIG. 5: The calibration of the extracted mtop value as a function of 
generated mtop is shown for the (a) e+jets and (b) /i+jets channels. 
The dependence is fitted to a linear function. Each point represents 
a set of 1000 pseudo-experiments for one of the fourteen (m, ,mf) 
combinations. Similarly, the pull widths, as defined in the text, are 
given for the (c) e+jets and (d) /i+jets channels. 
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FIG. 6: The normalized likelihood in Am'^^' after calibration via 
Eq. (15), together with a Gaussian fit, is shown for the (a) e+jets 
and (c) fl +jets events in data. The extracted Am'^"' values are indi- 
cated by arrows. The distributions in expected uncertainties af- 
ter calibration via Eq. (15) and correction for the pull width, obtained 
from ensemble studies using simulated MC events, is displayed for 
the (b) e+jets and (d) /x+jets channel. The observed S^fl values are 
indicated by arrows. 



C. Results 

With the calibration of Am and mtop, we proceed to extract 
Am and, as a cross check, mtop, from the data, as described 
in Sec. V. As indicated previously, the probabilities for the 
selected events are calculated using the ME method, and the 
likelihoods in Ain and mtop are constructed independently for 
the e+jets and /i+jets channels. 

The calibration of data involves a linear transformation of 
the uncalibrated axes of the likelihoods in Am and mtop to their 
corrected values, which we denote as Am'^^' and mf^p, accord- 
ing to: 



Ant 



cal 



Am-^0^" 



: Am 



„cal 
'top 



mtop -172.5 GeV-^, 



"Itop 





; '"top 



(15) 



172.5 GeV, (16) 



where the i^, are summarized in Table 4. The resulting like- 
lihoods for data, as a function of Am and mtop are shown in 
Figs. 6 and 7, respectively. 

After calibration, (Am) and (mtop) with their respective un- 
certainties and 5,„,^p, are extracted from the likelihoods as 
described in Sec. VI A. The uncertainties are scaled up by the 
average pull widths given in Table 4. The resulting distribu- 
tions in expected uncertainties are also shown in Fig. 6. 



The final measured results for Am and 



'«top 



are summarized 




FIG. 7: The normalized likelihood in mf^p after calibration via 
Eq. (16) together with a Gaussian fit for the (a) e+jets and (b) /i+jets 
channel. Arrows indicate the extracted mtop values. 



below according to channel, as well as combined: 



^ ^ ^ 1 Am 
e+jets, 2.6 fb"': 

'Wtop 



-1. 



^+jets, 2.6 fb 
£+jets, 2.6fb-' 



Am = 

'Wtop = 
Am = 



0.1 
173.9 
-0.5 
175.3 
-0.2 
174.7 



3.1 
1.6 
2.9 
1.3 
2.1 
1.0 



GeV 
GeV 
GeV 
GeV 
GeV 
GeV. 



(17) 



The uncertainties given thus far are purely statistical. The 
combined ^+jets results are obtained by using the canonical 
weighted average formulae assuming Gaussian uncertainties. 
We cross check the above values for mtop with those obtained 
from the absolute top quark mass analysis [43, 47] and find 
them to be consistent. 
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FIG. 8: Two-dimensional likelihood densities in m, and mf for the 
(a) e+jets and (b) /i+jets channels. The bin contents are propor- 
tional to the area of the boxes. The solid, dashed, and dash-dotted 
hnes represent the 1, 2, and 3 SD contours of two-dimensional Gaus- 
sian fits (corresponding to approximately 40%, 90% and 99% con- 
fidence level, respectively) to the distributions defined in Eq. (18), 
respectively. 



TABLE 5: Summary of systematic uncertainties on A/n. 



Source 



Uncertainty 
on Am (GeV) 



Modeling of detector: 




Jet energy scale 


0.15 


Remaining jet energy scale 


0.05 


Response to h and fight quarks 


0.09 


Response to h and b quarks 


0.23 


Response to c and c quarks 


0.11 


Jet identification efficiency 


0.03 


Jet energy resolution 


0.30 


Determination of lepton charge 


0.01 


ME method: 




Signal fraction 


0.04 


Background from multijet events 


0.04 


Calibration of the ME method 


0.18 


Total 


0.47 



As an additional cross check, we independently extract the 
masses of the top and antitop quarks from the same data sam- 
ple. The two-dimensional likelihood densities, as functions of 
nif and m,-, are displayed in Fig. 8. Also shown are contours 
of equal probability for two-dimensional Gaussian fits to the 
likelihood densities, where the Gaussian functions are of the 
form 



P{x,y) 



1 



r 1 1 

X exp 



- 21-p2 
2p{x-x){y-y)- 



OxOv 



(18) 



with x = mt and y = nif. The fits to data yield 



e+jets, 2.6 fb 



jU+jets, 2.6 fb 



-1. 



m, = 


173.8 


± 


1.5 


GeV 


m,- = 


173.8 


± 


2.0 


GeV 


P = 


-0.02 








nit = 


175.2 


± 


1.8 


GeV 


mf ~ 


175.5 


± 


1.5 


GeV 


P = 


-0.01. 









(19) 



The above uncertainties are again purely statistical; however, 
in contrast to Eq. (17), they are not corrected for pull widths in 
nit and nif. The correlation coefficients p are consistent with 
the absence of correlations. 

In Sec. VIII, we will combine the results for Am summa- 
rized in Eq. (17) with the previous measurement using 1 fb^' 
of integrated luminosity [12]. 



VII. SYSTEMATIC UNCERTAINTIES 

For the measurement of mtop we typically consider three 
main types of sources of systematic uncertainties [43]: 
(/) modeling of tt production and background processes. 



(n) modeling of detector response, and {Hi) limitations inher- 
ent in the measurement method. However, in the context of a 
Am measurement, many systematic uncertainties are reduced 
because of correlations between the measured properties of 
top and antitop quarks, such as, the uncertainty from the ab- 
solute JES calibration. Given the small value of the upper 
limit of (9(5%) already observed for |Am|/mtop, several other 
sources of systematic uncertainties relevant in the measure- 
ment of ffitop, such as modeling of hadronization, are not ex- 
pected to contribute to Am because they would affect t and t 
in a similar manner Following [48], we check for any effects 
on Am that might arise from sources in the latter category in 
Sec. VII C, and find them consistent with having no signif- 
icant impact. We therefore do not consider them further in 
the context of this measurement. On the other hand, we esti- 
mate systematic uncertainties from additional sources which 
are not considered in the m^p measurement, for example from 
the asymmetry in calorimeter response to b and b quark jets. 

Typically, to propagate a systematic uncertainty on some 
parameter to the final result, that parameter is changed in the 
simulation used to calibrate the ME method, and the Am re- 
sult is re-derived. If the change in a parameter can be taken 
into account through a reweighting of events, a new calibra- 
tion is determined using those weights and applied directly 
to data. When this procedure is not possible, a re-evaluation 
of event probabilities is performed for one sample of ti MC 
events coiTesponding to a particular choice of m, and nij clos- 
est to the most likely value according to our measurement, i.e. 
m, = nif =175 GeV, or, when no such sample of MC events 
with a changed parameter is available, m, = m,- = 172.5 GeV. 
Consequently, the results of ensemble studies are compared to 
those found for the default sample for the same values of nit 
and mf. 

The systematic uncertainties are described below and sum- 
marized in Table 5. The total systematic uncertainty is ob- 
tained by adding all contributions in quadrature. 
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A. Modeling of detector 

(i) Jet energy scale: As indicated in Sec. VI C, we use the 
absolute JES calibration of /tjEs = 1.018 ±0.008 deter- 
mined from data. To propagate this uncertainty to Am, 
we scale the jet energies in the selected data sample by 
fejEsilSD. 

(ii) Remaining jet energy scale: The systematic uncer- 
tainty on the absolute JES discussed above does not ac- 
count for possible effects from uncertainties on jet en- 
ergy corrections that depend on fijet and T7jet. To esti- 
mate this effect on Am, we rescale the energies of jets 
in the default // MC sample by a differential scale fac- 
tor ^(EjetjTjjet) that is a function of the JES uncertain- 
ties, but conserves the magnitude of the absolute JES 
correction. 

(iii) Response to b and light quarks: The difference in 
the hadronic/electromagnetic response of the calorime- 
ter leads to differences in the response to b and light 
quarks between data and simulation. This uncertainty 
is evaluated by re-scaUng the energies of jets matched 
to b quarks in the default tt MC sample. 

(iv) Response to b and b quarks: The measurement of Am 
can be affected by differences in the reconstruction of 
the transverse momenta of particles and antiparticles. 
A difference could in principle be caused by different 
Pt scales for /i+ and /i^. However, the data consist 
of an almost equal mix of events with opposite mag- 
net polarities, thereby minimizing such biases. We do 
not observe any difference in calorimeter response to e"*" 
and e^ . 

A systematic bias to Am can also be caused by dif- 
ferences in calorimeter response to quarks and anti- 
quarks. In the case of ti events, this bias could arise 
especially from a different response to b and ^-quarks. 
Several mechanisms could contribute to this, most no- 
tably a different content of K'^/K~ mesons, which have 
different interaction cross sections. In our evaluation of 
this systematic uncertainty, we assume that, although 
differences in response to b/b quarks are present in 
data, they are not modeled in MC events. We measure 
the difference of the calorimeter response to b quarks 
to that of b quarks, g = — using a "tag-and- 
probe" method in data. Namely, we select back-to-back 
dijet events, and enhance the bb content by requiring 
fo-tags for both jets. The tag jet is defined by the pres- 
ence of a muon within the jet cone, whose charge serves 
as an indication whether the probe jet is more likely 
to be a h or a /j-quarkjet. By evaluating the \pt\ im- 
balance between tag and probe jets for positively and 
negatively charged muon tags, we find an upper bound 
l-^hhl < 0.0042. Based on this result, we modify the 
default ti MC sample by re-scaling the momenta \p\ of 
b (fo)-quark jets by 1 T r %:t = 0.9979 (1.0021), and 
adjusting their 4-vectors accordingly. We repeat the en- 
semble studies after recalculating the probabilities for 



the modified sample and quote the difference relative to 
the default sample as a systematic uncertainty. 

(v) Response to c and c quarks: A difference in calorime- 
ter response to c and c quarks can potentially bias Am, 
since c quarks appear in decays of W+ bosons from t 
quark decays, and vice versa for c and /. It is exper- 
mentally difficult to isolate a sufficiently clean sample 
of cc dijet events, since it will suffer from considerable 
contributions from bb dijet events. However, the ma- 
jor underlying mechanisms that could cause a response 
assymetry, like, e.g., the different content of K^/K~ 
mesons, are the same, but of roughly opposite magni- 
tude between c and b quark jets, which would result in 
an anticorrelation. Based on the above, we assume the 
same upper bound \^cx\ < ^bh < 0.0042, and treat 

and g as uncorrelated. To propagate the sys- 
tematic uncertainty from ^c,c to we apply a simi- 
lar technique to that for the estimation of the systematic 
uncertainty due to different response to b and b quarks. 

(vi) Jet identification efficiency: DO uses scale factors to 
achieve data/MC agreement in jet identification effi- 
ciencies. To propagate to the Am measurement the ef- 
fect of uncertainties on these scale factors, we decrease 
the jet identification efficiencies in the default tt sample 
according to their uncertainties. 

(vii) Jet energy resolution: An additional smearing of jet 
energies derived by comparison of the pj balance in 
(Z ee) + Ijet events [49] is applied to all MC sam- 
ples in this analysis in order to achieve better data/MC 
agreement. To evaluate any effect from data/MC dis- 
agreement in jet energy resolutions on Am, we modify 
the default tt MC sample by varying the jet energy res- 
olution within its uncertainty. 

(viii) Determination oflepton charge: This analysis uses the 
charge of the lepton in ti candidate events to distinguish 
the top quark from the antitop quark. Incorrectly recon- 
structed lepton charges can result in a systematic shift 
in the measurement. The charge misidentification rate 
is found to be less than 1 % in studies of Z — > ee data 
events. To estimate the contribution of this uncertainty, 
we assume a charge misidentification rate of 1% for 
both e+jets and /J+jets final states and evaluate the ef- 
fects on Am resulting from a change in the mean values 
of the extracted m^^ and m?^. 

B. ME method 

(i) Signal fraction: The signal fractions / presented in Ta- 
ble 3 are changed by their respective uncertainties for 
each decay channel, and ensemble studies are repeated 
for all MC samples to re-derive the calibration for Am. 
The new calibrations are applied to data and the results 
compared with those obtained using the default calibra- 
tion. 



14 



(ii) Background from multijet events: In the calibration of 
this analysis, the background contribution to pseudo- 
experiments is formed using only W+jets events, as 
they are also assumed to model the small MJ back- 
ground from QCD processes and smaller contributions 
from other background processes present in the data. 
To estimate the systematic uncertainty from this as- 
sumption, we define a dedicated MJ-enriched sample 
of events from data. The calibration is re-derived with 
this background sample included in forming pseudo- 
experiments. 

(iii) Calibration of the ME method: The statistical uncer- 
tainties associated with the offset (i^o) and slope {^\) 
parameters that define the mass cahbration in Sec. VI B 
contribute to the uncertainty on Am. To quantify this, 
we calculate the uncertainty 5^,,, due to and 5^^ for 
each channel according to the error propagation formula 

and then combine the resulting uncertainties for the 
e-|-jets and /i-|-jets channels in quadrature. 

C. Additional checlis 



TABLE 6: Summary of additional checks for a possible bias on Am. 
None of those show any significant bias on Am. Note that the num- 
bers shown reflect an upper limit on a possible bias. This limitation 
is of statistical origin and due to the number of available simulated 

MC events. 



„ Change in Am 
^^ource ^^^^^ 


Modeling of physical processes: 




Higher-order corrections 


0.26 


ISR/FSR 


0.21 


Hadronization and underlying event 


0.23 


Color reconnection 


0.27 


^-fragmentation 


0.03 


PDF uncertainty 


0.10 


Multiple hadron interactions 


0.06 


Modeling of background 


0.07 


Heavy-flavor scale factor 


0.02 


Modeling of detector: 




Trigger selection 


0.07 


Zj-tagging efficiency 


0.25 


Momentum scale for e 


0.05 


Momentum scale for jU 


0.06 



(iii) Hadronization and underlying event: To check a possi- 
ble effect of Am from the underlying event as well as the 
hadronization models, we compare samples hadronized 
using PYTHIA with those hadronized using Herwig. 



We check for effects on Am from sources of systematic un- 
certainties considered in the mtop measurement [43] which are 
not expected to contribute any bias in the context of the mea- 
surement of Aw. For this, we follow the same approach as 
outUned at the beginning of this Section. We find the results 
of our checks to be indeed consistent with no bias on Am. 

The additional checks are described below and summarized 
in Table 6. Note that the numbers quoted merely reflect an 
upper bound on a possible bias, rather than any true effect. 
This limitation is statistical in nature and due to the number of 
available simulated MC events. Furthermore, if the difference 
between the central result and the one obtained for a check is 
smaller than the statistical uncertainty on this difference, we 
quote the latter. 



/. Modeling of physical processes 

(i) Higher-order corrections: To check the effect of 
higher-order corrections on Am, we perform ensemble 
studies using tt events generated with (/) the NLO MC 
generator mc@nlo [50], and {ii) the LO MC genera- 
tor ALPGEN, with Herwig [51] for hadronization and 
shower evolution. 

(ii) Initial and final- state radiation: The modeling of extra 
jets from ISR/FSR is checked by comparing pythia 
samples with modified input parameters, such as the 
±1SD changes, found in a study of DreU-Yan pro- 
cesses [52]. 



(iv) Color reconnection: The default PYTHIA tune used at 
DO (tune A), does not include explicit color reconnec- 
tion. For our check, we quantify the difference between 
Am values found in ensemble studies for ti MC sam- 
ples generated using tunes Apro and ACRpro, where 
the latter includes an explicit model of color reconnec- 
tion [53, 54]. 

(v) b-fragmentation: Uncertainties in the simulation of 
^-quark fragmentation can affect the measurement of 
mtop in several phases of the analysis, such as in b- 
tagging and in the /j-quark transfer functions used in the 
ME calculations. Such effects are studied in the context 
of Am by reweighting the simulated tt events used in 
the calibration of the method from the default Bowler 
scheme [55], which is fiined to LEP (ALEPH, OPAL, 
and DELPHI) data, to a tune that accounts for differ- 
ences between SLD and LEP data [56]. 

(vi) Uncertainty on PDF: The CTEQ6M [33] PDFs pro- 
vide a set of possible excursions in parameters from 
their central values. To check the effect on Am from 
PDFs, we change the default tt MC sample (generated 
using CTEQ6L1) by reweighting it to CTEQ6M, repeat 
the ensemble studies for each of the parameter varia- 
tions, and evaluate the uncertainty using the prescribed 
formula [33]: 

1 

Wdf = i|l?=2i[Am(5+)-Am(5r)]2|', 
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where the sum runs over PDF uncertainties for positive 

(5^) and negative (5^^) excursions. 

(vii) Multiple hadron interactions: When caUbrating the 
ME method, we reweight the luminosity profiles of our 
MC samples to the instantaneous luminosity profile for 
that data-taking period. For our check, we re-derive the 
calibration ignoring luminosity-dependent weights. 

(viii) Modeling of background: We check the effect of in- 
adequate modeling of background processes on our 
Am measurement by identifying distributions in the 
background-dominated £-|-3jets events that display 
only limited agreement between data and predictions 
from the sum of our signal and background models, as 
determined through a Kolmogorov-Smirnov test [57]. 
The calibration of the method is then re-done using 
W-l-jets events that are reweighted to bring the iden- 
tified distributions of predicted signal and background 
events into better agreement with data. 

(ix) Heavy-flavor scale-factor: As discussed in Sec. IV, a 
heavy-flavor scale-factor of 1 .47 ± 0.22 is applied to the 
W+bb+jets and W+cc +jets production cross sections 
to increase the heavy-flavor content in the alpgen 
W-njets MC samples. Moreoever, a scale factor of 
1.27 ±0.15 for the W+c+jets production cross section 
is obtained using mcfm. We re-derive the calibration 
with the heavy-flavor scale-factor changed by ±30% to 
check the magnitude of the effect on Am. 



2. Modeling of detector 

(i) Trigger selection: To check the magnitude the ef- 
fect from differential trigger efficiencies on Am, we 
re-derive a new Am calibration ignoring the trigger 
weights. 

(ii) b-tagging efficiency: We check the possibility of a bias 
in our Am measurement from discrepancies in the b- 
tagging efficiency between data and MC events by us- 
ing absolute uncertainties on the ^-tagging efficiencies, 
and account independently for possible discrepancies 
that are differential in rj and pr of the jet by reweighting 
the fc-tagging rate in simulated tt MC events to that ob- 
served in data. The total magnitude of a possible effect 
is determined by combining in quadrature excursions of 
Am values obtained with the modified calibrations for 
both absolute and differential changes. 

(iii) Momentum scale for electrons: DO calibrates the en- 
ergy of electrons based on studies of the Z ^ ee mass 
for data and MC events. We rescale the electron en- 
ergies in the default signal MC sample according to the 
uncertainties on the electron energy calibration to check 
the magnitude of the effect in the context of Am. 

(iv) Momentum scale for muons: The absolute momen- 
tum scale for muons is obtained from J/xj/ ^ jXjx and 



Z ^ li^i data. However, both linear and quadratic in- 
terpolation between these two points can be employed 
for the calibration. We check the effect of each extrap- 
olation on Am by applying the respective corrections to 
simulated tt MC events in the default sample, and find 
a larger shift in Am for the linear parametrization. 

VIII. COMBINING THE 2.6 fb" ^ AND 1 fb" ^ ANALYSES 

We use the BLUE method [58, 59] to combine our new 
measurement (Eq. 17) with the result of the analysis per- 
formed on data corresponding to Ifb"' [12]. The BLUE 
method assumes Gaussian uncertainties and accounts for cor- 
relations among measurements. 

For reference, we summarize the results obtained for 1 fb^ ' : 

e±jets, Ifb"': Am = 0.3 ± 5.0 (stat) GeV, 
;U±jets, Ifb"' : Am = 6.7 ± 4.7 (stat) GeV, 
^+jets, Ifb-i; Am = 3.8 ± 3.4 (stat) GeV. 

The 1 fb~^ analysis used a data-driven method to estimate 
systematic uncertainties from modeling of signal processes. 
This method did not distinguish between different sources 
of systematic uncertainties such as: (i) higher-order correc- 
tions, {ii) initial and final state radiation, [Hi) hadronization 
and the underlying event, and [iv) color reconnection. The 
above sources are studied in the context of the mtop measure- 
ment [43], but are not expected to contribute any bias to the 
measurement of Am. We cross-check their impact on Am in 
Sec. VllC, and find them consistent with no bias. Based on 
our findings, we do not consider any systematic uncertainties 
from modeling of signal and background processes. 

Two sources of systematic uncertainties from modeling of 
detector peformance (Table 5) are taken to be uncorrelated 
between the two measurements: JES and remaining JES. The 
rest are taken to be fully correlated. 

In the 1 fb~' analysis, a systematic uncertainty of 0.4 GeV 
from the difference in calorimeter response to b and b quarks 
was estimated using MC studies and checks in data. This 
systematic uncertainty has been re-evaluated using an entirely 
data-driven approach (item (iv) in Sec. VII A), and we there- 
fore use this new result for the analysis of the 1 fb" ' data. 
Furthermore, we now evaluated a systematic uncertainty from 
the difference in calorimeter response to c and c quarks, and 
propagate our findings to the 1 fb"^ analysis. 

AU other systematic uncertainties not explicitly mentioned 
above are taken as uncorrelated. 

The combined result for Am corresponding to 3.6 fb~' of 
data is 

Am = 0.84 ±1.81 (stat.) ±0.48 (syst.) GeV. (20) 

In this combination, BLUE determines a relative weight of 
72.8% (27.2%) for the 2.6 fb"' (Ifb^^) measurement. The 
X^/^DOF of the combination is 0.96. The combined likeli- 
hood densities for the two analyses are presented in Fig. 9 as 
functions of W; and Wf, separately for the e+jets and /i+jets 
channels. 
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FIG. 9: Combined likelihoods of the 2.6 fb"^ and Ifb"' measure- 
ments as functions of m, and mf in data for the (a) e+jets and 
(b) /j+jets channel. The bin contents are proportional to the area 
of the boxes. The solid, dashed, and dash-dotted lines represent the 
1, 2, and 3 SD contours of two-dimensional Gaussian fits defined in 
Eq. (18) (corresponding to approximately 40%, 90% and 99% con- 
fidence level, respectively) to the distributions, respectively. No pull 
corrections have been applied, and therefore the figures are for illus- 
trative purposes only. 



I. APPENDIX: GENERATION OF tt EVENTS WITH 



We briefly describe below the modifications to the stan- 
dard PYTHIA [32] code which were necessary to generate tt 
events with m, ^ mj. A new entry in the KF particle table is 
created for the t quark. The PYINPR subroutine is modified 
for use cases in which one of the tt production subprocesses 
(ISUB = 81,82,84,85) is called. The t quark is assigned as 
the second final-state particle whenever a t quark is selected 
as the first final-state particle. Furthermore, the ordering of the 
first and second final-state particles are swapped, as needed, 
in the subroutine PYSCAT. Additional changes are made in the 
subroutines PYMAXI, PYRAND, and PYRESD to set the lower 
limit on the combined masses of the W+ (W^) boson and b 
ib) quark to the t (f) quark mass. Finally, the subroutine PY- 
WIDT is modified to adjust the resonance widths and Ff as 
functions of and nif. 
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LX. CONCLUSION 
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about 3.6 fb"'. We find 
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ssl% of the mass of the top quark. 
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